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ABSTRACT 

Several recent numerical simulations of off-center cluster mergers predict that significant angular mo- 
mentum with associated velocities of a few xlO'^ km can be imparted to the resulting cluster. Such 
gas bulk velocities can be detected by the Doppler shift of X-ray spectral lines with ASCA spectrome- 
ters. Using two ASCA observations of the Centaurus cluster (Abell 3526), we produced a velocity map 
for the gas in the cluster's central regions. We also detected radial and azimuthal gradients in tem- 
perature and metal abundance distributions, which seem to be associated with the infalling sub-group 
centered at NGC 4709 (Cen 45). More importantly, we found a significant (>99.8% confidence level) 
velocity gradient along a line near-perpendicular to the direction of the incoming sub-group and with a 
maximum velocity difference of ^ 3.4±1.1 xlO^ km s~^. It is unlikely (P < 0.002) that the observed 
velocity gradient is generated by gain fiuctuations across the detectors. While the observed azimuthal 
temperature and abundance variations can be attributed to the interaction with Cen 45, we argue that 
the intracluster gas velocity gradient is more likely due to a previous off-center merging event in the 
main body of the Centaurus cluster. 

Subject headings: galaxies: clusters: individual (Abell 3526) — intergalactic medium — cooling fiows 
— X-rays: galaxies — 

ApJ in press 



1. INTRODUCTION 

Most of the visible mass in clusters is in the form of a 
very hot (T ~ 10^~* K) diffuse X-ray emitting plasma en- 
riched with heavy elements that permeates the cluster, the 
so-called intracluster medium (ICM). Galaxy clusters are 
believed to be formed by the infall/merging of smaller scale 
systems (bottom- up hierarchical scenario). The merging 
of sub-clumps induces several long lasting intracluster gas 
signatures that can be observed in the X-ray range, such 
as temperature and density inhomogeneities, destruction 
of cooling fiows and metal abundance gradients and gen- 
eration of gas bulk velocities. Spatial distributions of sur- 
face brightness and gas temperature can be directly com- 
pared to numerical hydrodynamical-f N body simulations 
of cluster mergers to provide a picture of the evolutionary 
stage of clusters (e.g. Evrard 1990; Katz & White 1993; 
Roettiger, Burns & Loken 1993,1996; Schindler & MuUer 
1993; Pearce, Thomas & Couchman 1994; Navarro, Frenk 
& White 1995; Evrard, Metzler, & Navarro 1996; Roet- 
tiger, Loken & Burns 1997, Ricker 1998; Takizawa & Mi- 
neshige 1998; Burns et al. 1999; Takizawa 1999, 2000 and 
references therein). 

ASCA observations have shown significant gas tempera- 
ture variations in many clusters, which are often attributed 
to recent mergers and this suggests that galaxy clusters are 
still forming. However, when determining the precise evo- 
lutionary stage of the intracluster gas, the distribution of 
gas temperature and surface brightness alone have limi- 
tations. This is because they represent only a 2-D pro- 
jection of a 3-D system. In order to reduce degeneracies 
and to determine the cluster evolution stage more precisely 
it is extremely important to obtain as much information 
as possible about the dynamics of the intracluster gas. 



Several recent numerical simulations of off-center cluster 
mergers predict that significant angular momentum can be 
imparted to the resulting cluster with velocities approach- 
ing a few thousand km s^^. Simulations also show that 
gas bulk fiows may last for several Gyr after the merger 
(e.g. Roettiger, Loken & Burns 1997, Ricker 1998; Roet- 
tiger, Stone, & Mushotzky 1998; Takizawa & Mineshige 
1998; Burns et al. 1999; Takizawa 1999, 2000; Roettiger 
& Flores 2000). Therefore, measurements of radial velocity 
distribution of the intracluster gas provides a very impor- 
tant additional source of information that, coupled with 
temperature and surface brightness distributions, can be 
directly compared to numerical simulations to determine 
the evolutionary stage of galaxy clusters. 

Radial velocity measurements of the intracluster gas can 
be obtained from X-ray observations, but require accurate 
determinations of spectral line centroids. The precision 
with which a line centroid can be measured is a func- 
tion of the square root of the number of line photons, the 
FWHM of the instrument, the line energy and, in practice, 
it strongly depends on knowledge of the instrumental gain 
(conversion between pulse height and energy). The instru- 
ments on-board ASCA have the capability of measuring 
directly such ICM bulk velocities. Furthermore, ASCA's 
spectrometer's gain calibration is already well-understood 
as to allow reliable velocity measurements. 

Recently, using the GISs on-board ASCA, we have de- 
tected a velocity gradient in the intracluster gas of the 
Perseus cluster. The spectral analysis of 8 different ASCA 
pointings encompassing a region of ~ 1 h^Q^ Mpc radius 
around the center of the Perseus cluster (Dupke & Breg- 
man 2000, 2001), showed a significant velocity gradient 
of > 1000 km s"^Mpc-i at the 90% confidence. How- 
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ever, the GISs have worse energy resolution than the SISs, 
which severely limited the ability to precisely constraint 
the velocity gradient amplitude. Ideally, one would like 
to use the higher energy resolution of the SISs to measure 
ICM bulk velocities. 

The Centaurus cluster provides us with an excellent op- 
portunity to measure intracluster gas velocities more pre- 
cisely. It is a bright, nearby cluster, with high metal abun- 
dances (well isolated spectral lines) and shows signs of in- 
teraction with an infalling subgroup (Ccn 45) located ^ 
13-18' from the center, which is associated with the galaxy 
NGC 4709. Furthermore, it has been an ASCA target for 
a long exposure 90 ksec) taken in 1-CCD mode with 
the best calibrated CCDs. These observational character- 
istics make Centaurus an optimal target for gas velocity 
measurements. 

In this paper we analyze the central regions of Centaurus 
and determine the velocity structure of this system using 
data from two different ASCA pointings. In sections 2 & 
3 we describe the general characteristics of the Centaurus 
cluster, of the ASCA pointings analyzed and of data re- 
duction procedures. In section 4 we investigate whether 
the velocity differences between Cen 30 and 45 previously 
observed optically by other authors can be corroborated 
by X-ray velocity analysis using the GISs. In section 5 
we present the spectral fitting results of SIS data for the 
observation with the longest exposure (taken in 1995) and 
estimate the effects of gain variations in our results as a 
function of frequency and position. In section 6 we show 
the results of the spectral analysis of another ASCA point- 
ing taken 2 years earlier, with a shorter exposure. In sec- 
tion 7 we briefly discuss the nature of the small and large 
scale bulk motions detected in Centaurus. 

2. THE CENTAURUS CLUSTER 

Centaurus (Abell 3526) is classified as Bautz-Morgan 

type I, it has an optical redshift of 0.0104 and is one of 
the closest X-ray bright clusters of galaxies. It has a small 
cooling flow and the estimated accretion mass rate is < 30- 
50 solar masses per year (e.g. Matilsky, Jones & Forman 
1985; Allen & Fabian 1994; White, Jones & Forman 1997; 
Peres et al. 1998). Its X-ray emission is relatively smooth, 
slightly elliptical in shape and is strongly peaked on the 
cD galaxy NGC4696. The average temperature of the gas 
in the cluster is ~ 3.5 keV (e.g. Mitchell & Mushotzky 
1980; Matilsky et al. 1985; Thomas, Fabian fc Nulsen 
1987: Edge & Steward 1991). ASCA CIS & SIS analyses 
of the central region of Centaurus have shown evidence of 
a strong central metal abundance enhancement (Allen & 
Fabian 1994; Fukazawa et al. 1994) varying from superso- 
lar in the very central regions down to 0.3 solar at ~ 13' 
(1' at Centaurus distance is equivalent to ^ 19 h^o'^kpc). 
Furthermore, the region where metal abundances are en- 
hanced is dominated by SN la ejecta (Allen et al. 2000) 
indicating the prc^scmcc^ of a radial "chemical gradient" . 

Early studies of galaxy velocities in the Centaurus clus- 
ter have shown the existence of bimodality in velocity 
space (Lucey, Currie & Dickens 1986a,b). Two galaxy 
groups are clearly separated: The main group (Cen 30) 
which is centered on the cD galaxy (NGC 4696) shows 
an average radial velocity of '^3000 km s~^ and a veloc- 

^ http://heascirc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 



ity dispersion of ~590 km s~^. The second group (Cen 
45) is associated with the galaxy NGC 4709 at ^15' from 
NGC 4696. It has an average radial velocity of ~ 4600 
km s~^ and a velocity dispersion of ~ 280 km s~^. Lucey 
et al. (1986a, b) concluded that the Cen 45 and Cen 30 are 
at the same distances and that Cen 45 is being accreted 
by the Centaurus cluster, based on a detailed analysis of 
the color-magnitude, luminosity function and galaxy ra- 
dius distribution. Improved redshift measurements of the 
central regions of Centaurus confirmed the galaxy velocity 
bimodal distribution: the "true" cluster Cen 30 having an 
average velocity of 3397±139 km s~^ and a velocity dis- 
persion of 933±118 km s-\ and Cen 45 having an average 
velocity of 4746±43 km s~^ and a velocity dispersion of 
131±43 km s-i (Stein et al. 1997). 

The X-ray surface brightness contours of Centaurus are 
slightly elliptical and NGC 4709 is located near the major 
axis (Allen & Fabian 1994), suggesting that the elongation 
towards Cen 45 can be associated with the elongation of 
the gravitational well along the collision axis. More re- 
cently, Churazov et al. (1999) has derived an X-ray tem- 
perature map of the Centaurus cluster using the ASCA 
satellite and found a region where the intracluster gas is 
significantly hotter (^5 keV) than the surroundings (^^3.5- 
4.0 keV). The region with higher gas temperature is appar- 
ently associated with Cen 45 and indicates the presence of 
shocks due to the infall of subcluster Ccn45. 

3. DATA REDUCTION & ANALYSIS 

ASCA carries four large-area X-ray telescopes, each with 
its own detector: two Gas Imaging Spectrometers (GIS) 
and two Solid-State Imaging Spectrometers (SIS). Each 
GIS has a 50' diameter circular field of view and covers a 
energy range of 0.8 10 keV, while each SIS has a 22' square 
field of view and covers an energy range of 0.4-10 keV. For 
all pointings analyzed in this work we selected data taken 
with high and medium bit rates, with cosmic ray rigidity 
values > 6 GeV/c, with elevation angles from the bright 
Earth of > 20°, and from the Earth's limb of > 5° (GIS) 
or > 10° (SIS); we also excluded times when the satellite 
was affected by the South Atlantic Anomaly. Rise time re- 
jection of particle events was performed on GIS data and 
hot and flickering pixels were removed from SIS data. We 
estimated the background from blank sky files provided by 
the ASCA Guest Observer Facility. In the spectral fittings 
we used XSPEC vlO.O software (Arnaud 1996) to fit the GIS 
& SIS spectra separately and simultaneously. The spectra 
were fitted using the mekal and Vmekal thermal emission 
models, which are based on the emissivity calculations of 
Mewe & Kaastra (Mewe, Gronenschild & van den Oord 
1985; Mewe, Lemen & van den Oord 1986; Kaastra 1992), 
with Fe L calculations by Licdahl et al. (1995). Abun- 
dances are measured relative to the solar photospheric val- 
ues (Anders & Grevesse 1989), in which Fe/H=4.68 x 10"^ 
by number. Galactic photoelectric absorption was incor- 
porated using the wabs model (Morrison & McCammon 
1983); Throughout the paper the hydrogen column den- 
sity is assumed to be fixed to the Galactic value, i.e., 
8.07xl020cm-2 (Dickey & Lockman 1990 - W3nH tool^), 
unless stated otherwise. Spectral channels were grouped 
to have at least 25 counts/channel. Since 1994, there have 
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been increasing discrepancies between SIS and GIS in the 
low energy end of the spectra , due to a decrease in X-ray 
efficiency in the SISs^; Therefore, energy ranges for both 
instruments were restricted to 0.8-10 keV, in all spectral 
fittings. 

The central regions of the Centaurus cluster were ob- 
served by ASCA in June of 1993 (P93) and on July of 
1995 (P95) for ~ 24 ksec and 87 ksec, respectively; the 
characteristics of these two pointings are shown in Table 
1. The original energy resolution of the spectrometers on- 
board ASCA were 8% and 2% at 5.9 keV for GISs and 
SISs, respectively. The SISs energy resolution has been 
steadily degrading with time (e.g. Dotani et al. 1997) due 
to the increases in dark current and decline in X-ray sig- 
nal charge transfer efficiency (Rasmussen et al. 1998). For 
the best pointing analyzed in this work (P95) the energy 
resolution is ~ 2.9% and 3.1% at 6.0 keV for SIS & 1, 
respectively. Furthermore, for P95 the SIS observations 
were taken in 1-CCD mode using the best calibrated chips 
of each spectrometer (SOCl and S1C3). 

Since there is a mild cooling flow at the center of Cen- 
taurus we added a cooling flow component (Mushotzky & 
Szymkowiak 1988) to the mekal thermal emission model 
in the central region, to compare the temperature of the 
hot component in the central region with that of the outer 
regions. We tied the maximum temperature of the cooling 
flow to the temperature of the thermal component, and 
we fixed the minimum temperature at 0.1 keV. The abun- 
dances of the two spectral components (mekal and cf low) 
were tied together. We also applied a single (but variable) 
absorption to both spectral components. 

We have noticed an unusual gain variation in the SISs 
when comparing spectral fittings including the FeK com- 
plex to those excluding it. It is probably due to the as- 
sumption that the charge transfer inefficiency is indepen- 
dent of X-ray energy during the initial gain standard cor- 
rection procedure (see below). Typically, spectral fittings 
that excluded the FeK line complex showed higher best-fit 
redshift values. Since we are looking for relative velocity 
measurements this effect should not affect our analysis. 
However, we also show the results from spectral fittings 
excluding the FeK complex, where the energy range con- 
sidered goes from 0.8-6.0 keV and are denoted by a sub- 
index FeL. 

For P95 we noticed significant discrepancies of best-fit 
redshifts for some regions (e.g. P8) between SIS & 1 
spectral fits. This is probably due to different spatial dis- 
tribution of gain fluctuations in different chips. Therefore, 
we avoided performing simultaneous spectral fittings us- 
ing these two spectrometers. The statistics are sufficiently 
good that radial velocity distributions can be measured 
precisely on each instrument separately. For P93 how- 
ever, the statistics are not good enough to notice signifi- 
cant changes of the best-fit parameters when spectra from 
SIS & 1 are fitted separately. Therefore, for that point- 
ing, we show only the results of joint spectral fittings. 

For the analysis of the cluster's central regions with the 
SISs we selected a set of 9 different extraction regions 
within 8' of the cluster's center using the SISs. The very 
central region included the cluster core (PO), and the other 
regions were distributed symmetrically around PO, cen- 

^ heascirc.gsfc.nasa.gov/docs/asca/watchout.html, also Hwang et al. 



tered 5' away and with a radius of 3'. These regions are il- 
lustrated in Figure 1, where they are overlaid on the X-ray 
surface brightness contours. This set of extraction regions 
was applied to both pointings (P95 & P93). We applied a 
correction for the position angles between these pointings 
(19.6°), so that each extraction region corresponded to the 
same position in the sky for both observations. 

Since P93 was taken in 4-CCD mode and the bulk of the 
cluster emission was extended over 2 CCDs in both SIS 
& 1, we have taken into account the overall interc;liip gain 
differences when comparing spectra from different CCDs 
(see below). In doing so, if an extraction region crossed 
over 2 CCDs we considered only the fraction of that region 
that included the largest number of counts in one CCD. 
This procedure minimizes the errors involved in the pro- 
cess of gain correction. Therefore, when comparing the 
results from pointings P93 to P95 one should notice that 
the ccntroids of the extraction regions may be slightly (< 
1') shifted with respect to one another for P93 regions: PI 
(shifted towards P2), P4 (shifted towards P3), P5 (shifted 
towards P6) and P8 (shifted towards P7). 

For the GIS analysis of the outer (>8') parts of Cen- 
taurus, we selected two sets of extraction regions: in the 
first set, the region encompassing Gen 45 (circle centered 
11' from NGC 4696 towards Cen 45 with a 8' radius) was 
compared to the average. The second set of regions was 
chosen analogously to that for the SIS analysis of the cen- 
tral 8' and it included a central region (0-5') and 8 large 
regions symmetrically distributed around the center, at a 
distance of 12' from the center and with a radius of 5'. 
Since the GISs have a poorer energy resolution than the 
SISs, we joined the 8 outer regions by adjacent pairs (1-2, 
3-4, 5-6, 7-8) in order to improve photon statistics. These 
regions are also illustrated in Figure 1. The objective of 
the analysis of the outer regions was to check for corrobo- 
ration in X-ray frequencies of the velocity difference seen 
in the optical analyses between Cen 45 and Gen 30. 

4. GIS GLOBAL VELOCITY MEASUREMENTS OF CEN 45 

As we mentioned previously, the galaxy velocity distri- 
bution in Centaurus is bimodal (Lucey, Currie & Dick- 
ens 1986a, b). One of the velocity components is associ- 
ated with the main body of the cluster (CenSO) and has 
a radial systemic velocity of ^ 3000 km s^^ and another 
component associated with the galaxy NGC 4709 ~ 15' 
away from Cen 30 has a radial systemic velocity of ~ 
4500 km s^^ and is generally denoted by Cen 45. We 
examined whether the optical velocity differences have an 
X-ray counterpart. We applied a similar technique that 
was used in the analysis of velocity measurements in the 
Perseus cluster with the GISs on-board ASCA (see Dupke 
& Bregman 2001 for details). 

The most obvious spatial configuration to test for veloc- 
ity differences between the region associated with Cen45 
and Cen 30 is to extract spectra from a circular region 
large enough to encompass the bulk of Cen 45 (11' away 
from the center with a 8' radius - denoted C45) and com- 
pare the best-fit redshift values to the average (total field- 
of-view minus C45). The best-fit parameters of the spec- 
tral fittings of the C45 region and of the average region 
are shown in Table 2. C45 region is significantly hotter 
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(T=5.47±0.27 kcV) than the average gas temperature of 
Centaurus (3.44±0.03 keV), and has a lower metal abun- 
dance (0.5±0.09 Solar) than the average (0.77±0.03 So- 
lar), consistent with a central negative radial abundance 
gradient. The best-fit redshift value for C45 is found to be 
3.27li:^^ X 10-^ which is significantly (>90% confidence 
limits) higher than the average redshift of the rest of the 
cluster 1.81+2;^'^ x IQ-^. 

The errors cited above are spectral fitting errors and 
do not include any uncertainties related to gain variations 
across the detectors. In order to estimate the latter we 
used a similar procedure to that used previously in the GIS 
velocity analysis of Perseus (Dupkc & Brcgman 2001). We 
used the GIS gain map determined from the Cu-K line at 
8.048 keV (Idesawa et al. 1995) to estimate the velocity 
differences expected purely from residual gain variations 
across the detector. We compared the average gain values 
of the C45 region for each instrument with that of the field- 
of-view excluding C45. The gain variations across each 
detector are corrected in-house with the ftool ASCALIN 
V.9.0t. If we assume that there are still significant gain 
fluctuations across the detector and, conservatively, that 
these fluctuations are on the same order as those of the 
pre-processed data then the expected difference between 
the two regions due purely to gain variations is ~1870 
kms""* (C45 having the higher radial velocity). Therefore, 
the "real" velocity difference (not directly associated with 
gain variations) is ^ 2.6±2.G x 10"^ km (l-cr errors), 
which is consistent with the optical velocity measurements 
of ~ 1500 km s-^ 

The gain scatter across the C45 region in GIS 2 & 3 is 
relatively high (on the same order of the velocities that 
we are trying to measure) lowering the significance of the 
observed velocity differences. Therefore, we looked for an- 
other spatial extraction configuration that minimized the 
gain fluctuations inside the extraction regions but still en- 
compassing Cen 45 and with enough photon statistics to 
measure velocities correspondent to the optical measure- 
ments. The resulting configuration chosen is also shown 
in Figure 1 and it consists of 4 regions symmetrically dis- 
tributed aromid the central region (circle centered on the 
X-ray peak with a 5' radius). The four regions are pairs of 
adjacent circular regions (12' away with 5' radius). They 
are labeled as J12, J34, J56 & J78. This configuration 
shows a reduction of intrachip gain scatter by ~ 17% for 
GIS 2 and ^32% for GIS 3 with respect to the previous 
configuration. The high temperature region detected by 
Churazov et al. (1999) associated with Cen 45 is located 
mostly within ,156. The best-fitting parameter values for 
this configuration are also shown in Table 2. It can be 
seen that the only region which shows significantly differ- 
ent redshifts (>90% confidence) with respect to the central 
region is J56 an the corresponding velocity difference is ~ 
1.77±0.96 X 10-'' km s^^ (l-cr errors). 

The results above are consistent with the optical mea- 
surements and indicate that significant velocity gradients 
in clusters can be detected even with the GISs, which have 
worse energy resolution than the SISs on-board ASCA. In 
the next sections we use the SISs to analyze the very cen- 
tral regions of the Centaurus cluster. Besides having a 
better energy resolution (by a factor of ;^2.5), the SIS 
spatial gain fluctuations are more stable than that of the 



GIS. Therefore, radial velocities can be measured with a 
substantially higher precision. 

5. SIS ANALYSIS OF THE CENTRAL REGIONS OF 

CENTAURUS 

5.1. SIS Observation of July 1995 (P95) 

P95 is the best pointing analyzed in this work and prob- 
ably one of the best pointings in the ASCA archive for 
the analysis of gas bulk motions. The combination of 
long-exposure (Table 1), clock mode (1-CCD), and clus- 
ter's center positioning in the CCDs (roughly at the center 
of SOCl and S1C3) makes this pointing near-optimal for 
velocity studies with ASCA. The best-fit values for tem- 
perature, abimdance and redshift obtained from spectral 
fittings of the central regions are shown in Table 3 and 
plotted in Figures 2a,b as a function of the position angle 
for SIS & 1, respectively. The indicated errors in Figures 
2a, b are 90% confidence limits. 

It can be seen from Figures 2a, b that the temperature 
distribution is consistent with the presence of a cooling 
flow (with a best flt mass deposition rate <^ 10 solar masses 
per year). The best-fit temperatures are found to be typi- 
cally slightly higher using SIS 1 data than those obtained 
with SIS 0. The temperature in the center is ~2.68±0.05 
keV and shows a positive gradient outwards up to a max- 
imum value of ~3.59±0.09 keV and ~3.85±0.12 keV for 
SIS &1, respectively. The 90% confidence limits for the 
best-fit values in the central region are indicated by the 
two horizontal lines in Figures 2a,b. Furthermore, the ob- 
served temperature gradient is not purely radial. One can 
notice a significant azimuthal temperature gradient, which 
is roughly symmetric and with a maximum towards the di- 
rection of Cen 45, which is the region in between the two 
vertical lines in all plots. The azimuthal temperature dis- 
tribution has only one axis of symmetry and, therefore, 
is not consistent with the "cross" pattern shape of ASCA 
PSF scattering. 

The metal abundance distribution is consistent with the 
existence of a central abundance gradient (Figure 2a,b and 
Table 3). In general, if different elemental metal abun- 
dances were tied such as to produce solar abundance ra- 
tios the fits were significantly worse (often with ~ 1-3). 
There was a significant improvement in the spectral fit- 
tings when individual metal abundances were let free to 
vary (x^ w 1), resulting in strongly non-solar abundance 
ratios. This is consistent with the observed dominance of 
SN la ejecta in the central regions of this cluster (Allen et 
al. 2000). Furthermore, letting the individual elemental 
abundances free to vary did not cause significant changes 
in the best-fit values of temperature, Fe abundance or red- 
shifts. Therefore, we report here only the spectral fittings 
for an absorbed (Galactic) Vmekal (abundances free to 
vary) spectral model (with an extra cooling flow compo- 
nent for the central region). The metal abundance within 
the inner 3' was found to be 1.13±0.05 solar and 1.23±0.06 
solar, for the SIS & 1, respectively. The 90% confidence 
limits for the abundance in the central region are shown 
by the horizontal lines in Figures 2a, b. One can see also 
indications of an azimuthal metal abundance gradient. In 
the outer 5' the metal abundance drops significantly to- 
wards the direction of Cen 45 (P3, P4, P5, P6 & P7 for 
SIS and P4, P5 & P7 for SIS 1). The minimum metal 
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abundance values were reached around pointings P4 & P5 
in both detectors, where the best-fit value is ~0.86±0.1 
solar. 

The azimuthal distributions of redshifts are shown in 

Figures 2a, b (bottom plot) as a function of the position 
angle. The original SIS best-fit redshift values were too 
close to so that full determination of 90% confidence er- 
rors was compromised in XSPEC. Therefore, all redshifts in 
that chip were boosted up by a slope correction of 0.9%"^. 
This redshift-boosting does not affect our results since we 
are not interested in the absolute redshift values, but only 
in the redshift differences between regions. SIS 1 redshifts 
are displayed directly (without any gain boosting). 

From Figure 2a we can see that the SIS data shows a 
very significant velocity gradient. Regions P3 & P4 have 
significantly lower velocities than regions P6, P7 & P8 by 
~2.67±0.21 X 10^ km s"^ for SIS and 2.73±0.27 x 10^ 
km s""'^ for SIS 1 (1-(t errors). Both velocity-discrepant 
sets of regions have significantly different velocities (on 
the order of ^1350 km s~^ at > 90% confidence level) 
with respect to the clusters center. Although the velocity 
gradient is also seen in SIS 1 data, its significance is lower 
than that obtained using SIS data. This is due to the sig- 
nificantly smaller effective number of counts of SIS 1 with 
respect to SIS (30% less counts) and to the larger gain 
spatial variability of SIS 1. Figure 2b shows the azimuthal 
distribution of best-fit redshifts obtained using SIS 1 data. 
It can be seen that for the "low velocity" regions (P3 & 
P4) the redshifts are significantly different from that of 
the central region and also from most of the regions in 
the "high velocity" hemisphere (P6 & P8). There is one 
region (PS) for which analysis of SIS & 1 spectra gives 
discrepant results. Although the distribution of radial ve- 
locities determined using SIS 1 are less constrained than 
that determined using SIS 0, the overall azimuthal distri- 
butions of velocities are consistent with one another. The 
velocity gradient found with SIS 1 is maxim,um, along the 
same direction (P4-P6) as that determined from SIS and 
has the same magnitude 380 km s~^arcmin~^). The ve- 
locity differences are unlikely caused by spectral modeling 
effects since regions with discrepant velocities have similar 
temperatures and metal abundances. 

To assess the significance of the velocity differences de- 
scribed in the previous paragraph, we chose two regions 
where we consistently found the maximum velocity differ- 
ences. P4 and P6, and we applied the F-test. We simul- 
taneously fitted spectra from SIS & 1 for each of the 
two regions and analyzed variations associated to the 
change in the number of degrees of freedom. We com- 
pared the of fits which assumed the redshifts to be the 
same in the two projected spatial regions P4 & P6, i.e., 
zsopi = zslpi = zsope = zsipB, to that of fits which 
allowed the redshifts in the two regions to vary indepen- 
dently, reducing the number of degrees of freedom by one. 
For the latter we tied the redshifts within each of the two 
regions together, i.e., (zsop4 = zsip^) + (zso^-e = ^sxp^)-, 
where zsipj is the redshift of region j with instrument 

where the new hne energy is related to the old energy in the standard way, i.e., 

www. astro. isas.ac.jp/'~dotani/rdd. html 
^ heasaJTC. gsfc.nasa.gov/docs/asca/newsletters/sis_calibration5. html 
^ heasaj:c.gsfc.nasa.gov/docs/asca/newsletters/sis_performance2.html 
^ heasaj:c.gsfc.nasa.gov/docs/asca/newsletters/sis_calibration5.html 



SIS i. This is justified by the lack of significant changes 
in the best-fit redshifts obtained by different instruments 
for the same region. The difference between the of 
these two fits must follow a x^ distribution with one 
degree of freedom (e.g. Bevington 1969). The value 



, where n is the number of fitting 



terms, is related to the probability that the additional 
term significantly contributes to the reduction of the x^ 
statistic. Comparison of the F{n) from the fits with locked 
and unlocked redshifts to tabulated values of F indicates 
that the best-fit redshifts obtained for P4 & P6 are dis- 
crepant at > 99.99% level. We show in Figure 3 the 90% 
and 99% confidence contours for two interesting parame- 
ters (redshifts) of regions P4 and P6 as well as the line 
correspondent to equal redshifts. 

5.2. SIS Gain Variations & Redshift Dependence on 

Gain 

The SIS performance has gradually degraded due to the 
accumulated energetic particle radiation damage, which is 
manifested by an increase in residual dark-current distri- 
bution (RDD) and in signal charge transfer inefficiency 
(CTI). Although the RDD effect is neghgible in 1-CCD 
mode (Dotani 2000 ^) the increase of the CTI reduces the 
apparent X-ray energy and thus causes a systematic vari- 
ation of gain with time (Dotani et al.^, Rasmussen et al. 
^). This effect is different for different chips resulting in 
enhancements of interchip global gain variations. 

Correction for the above mentioned effect have been car- 
ried out using the Ni background fluorescence line origi- 
nated from the kovar (iron- nickel-cobalt alloy), which cov- 
ers the frame store region of the chips and also using ob- 
servations of Cas A. To investigate the time dependence 
of the average CTI, the center energy of the Ni line has 
been recorded every 3 months from standard chips (SOCl 
and S1C3) in 1 CCD mode (Dotani et al.''). The decrease 
of the line energy is well approximated by a linear func- 
tion of time. Furthermore, the gain variation with time 
was observed to vary non-uniformly for different spectral 
lines after CTI correction. This may be a consequence of 
the assumption that the CTI is independent of X-ray en- 
ergy as well as from scatter of the line energies among the 
pointing positions. 

The event files analyzed in this work were corrected by 
the tool ASCALIN V0.9t, using the most recent gain cal- 
ibration file sisph2pi_110397.fits. The line center energies 
are different from chip to chip by at most 2% (interchip 
gain variation). Siiic;e we are showing results for SIS & 
1 separately, we do not worry about global gain variation 
differences between these detectors. However, we need to 
estimate the gain variation across each CCD (intrachip 
gain variation). We did so by analyzing the gain scatter 
across the detectors measured using Cas A data. After the 
standard correction analysis is completed, the systematic 
errors on the FeKa line centers are reduced to ~ 0.29% 
(870km/s) for SOCl and ~ 0.28% (840 km/s) for S1C3. 

- intercept 



slope 
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The values mentioned above are based on a more recent 
preliminary calibration analysis performed by Dotani and 
kindly provided to us by K. Mukai ^. Using the older cali- 
bration references cited in the previous paragraphs^ would 
give a a gam of 0.22% and 0.38% for SIS & SIS 1, respec- 
tively. This discrepancy does not change any of the results 
obtained in this paper. However, the larger observed un- 
certainties associated with the determination of best-fit 
redshifts in SIS 1 suggests that again for SIS 1 should be 
greater than that for SISO. 

We tested the sensitivity of our observations to possible 
residual gain variations across the SISs using Monte Carlo 
simulations. Supposing the rcdshift to be constant for the 
two regions where wc find maximum velocity differences 
(P4 & P6), we generated fake spectra for both spectrom- 
eters and compared the best-fit redshift differences. Then 
we calculated the probability to find the same redshift dif- 
ferences (or greater than) that wc observed in the real 
regions , for SIS & 1. The procedure used is completely 
analogous to that used in testing the effects of residual 
gain variations across the GISs in the velocity measure- 
ments for the Perseus cluster (Dupke & Bregman 2000, 
2001) and wc refer the reader to that paper for details. 

We simulated 1000 SIS & 1 spectra corresponding to 
the real observation using the same spectral model, pa- 
rameters and responses correspondent to the real regions. 
The simulated spectra were then grouped, fitted and had a 
gain uncertainty added to their best-fit redshift values. We 
assumed that the gain variations follow a Gaussian distri- 
bution with a standard deviation {a gain) , which is different 
for each spectrometer, and a zero mean. We adopted as 
our 1-(T gain variations [ugain) f^or the SIS & 1 the val- 
ues described above, i.e., 0.29% and 0.28%, respectively. 
We then calculated the probability of finding the same (or 
greater) redshift differences between the two simulated re- 
gions that wc found in the real ones, i.e., Az = 0.009; that 
probabihty is found to be <^ 0.002. 

In Figure 4 we plot the dependence of the chance prob- 
ability (described in the previous paragraph) on the as- 
sumed again of the spectrometers (for the purpose of illus- 
tration again for SIS is assumed to be the same as that for 
SIS 1, which roughly represents the real case). In Figure 5 
we also indicate the value for again estimated in the pre- 
vious paragraphs (~0.28%). In our real case we observed 
not only two but at least three regions (more than three 
if SIS data alone is considered) with discrepant redshifts 
(P3, P4 & P6). Therefore, the significance of the veloc- 
ity gradient found in Centaurus obtained by the method 
described above (99.8%) should be taken as a lower limit. 

5.3. Spectral Fittings Without the FeK Line Complex 

The velocity measurements presented in the previous 
section were derived from spectral fittings of SIS & 1 
encompassing the usable energy range of the spectrome- 
ters, i.e., from 0.8 to 10 keV, and arc mainly driven by the 
FcKq line complex. Since the Centaurus cluster has very 
high (^solar) metal abundances it allows for the possibility 
of reliably measuring velocities using other spectral lines, 
or at least, excluding the FeK line complex energy range. 
When we fitted the spectra from SIS & 1 for different 

^ see www.astro.Isa.umich.edu/~rdupke/dotani.ps 

^ heasaj:c.gsfc.nasa.gov/docs/asca/newsletters/rad_dam_sis3.html 



energy ranges we noticed a systematic gain discrepancy 
between spectral fittings that included FeK line complex 
from those excluding it. Excluding energies corresponding 
to both FeK and FeL complexes in the spectral fittings pro- 
duced best-fit velocities (determined from spectral lines of 
other individual elements) that were more consistent with 
those fits that included the FeL line complex than with 
those that included the FeK line complex. When the FeK 
lines were included the measured redshifts were consis- 
tently lower. As mentioned previously, this effect is most 
likely due to the assumption that the CTI is independent 
of X-ray energy during the in-house calibration process. 

In order to check how the observed velocity distribu- 
tions were affected by the effect mentioned in the previous 
paragraph we performed a series of redshift measurements 
excluding the FeK region (<6 keV). The spectral fitting 
characteristics were completely analogous to that used in 
the last section. The resulting best-fit values for temper- 
ature, metal abundance and radial velocities are listed in 
Table 4 and the redshift distributions for SIS & 1 are 
shown in Figures 5a, b. The best-fit temperatures found 
when the FeK line complex is excluded from the analy- 
sis are, in general, slightly higher than, but similar to, 
those obtained using the full useable energy range. The 
azimuthal temperature gradient is still marginally seen. 
The measured temperatures in the center are 2.84±0.03 
keV (SISO) & 2.771^ °? keV (SISl), rising to 3.70±0.12 
kcV (SISO) & 4.11±0.16 keV (SISO) towards the direction 
of Cen 45. The abundance measurements are not well con- 
strained enough to detect a central metal enhancement and 
are consistent with a flat proflle (~ 1.25±0.25 Solar) with 
a large dispersion. 

It can be seen from Figures 5a, b that, although the best- 
fit redshifts are typically higher than those in Figures 2a,b, 
the behavior of the azimuthal distributions with and with- 
out the FeK complex energy range are very similar. A sig- 
nificant (> 90% confidence) velocity gradient is also found. 
The velocity distribution reaches a minimum around re- 
gions P3 & P4 and a maximum at region P6 & P7. The 
velocity gradient is found to be steeper for SIS 0. The cor- 
responding maximum velocity differences are ^ 6.3±1.2 x 
10^ km s-i for SIS and ~'3.18±0.96 x 10^ km for 
SIS 1 (1-cr errors). 

5.4. SIS Observation of June 1993 (P93) 

In the previous sections we have shown that the velocity 
gradient in the ICM of the Centaurus cluster is detected 
with high significance using SIS & 1 data including and 
excluding the FeK line complex. We have also shown that 
it is not due to gain variations within the CCDs. Since 
ASCA has observed the Centaurus cluster several times, 
we can test the reliability of the velocity measurements by 
performing the same analysis for another pointing, with 
different observing characteristics (e.g. observing date, 
CCD clock mode, position on the CCD, etc). Out of the 
two archived pointings in the ASCA archive that cover the 
whole central regions of Centaurus, we chose the one with 
larger spatial coverage and exposure time. We refer to this 
pointing as P93 and list its characteristics in Table 1. 
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P93, however, has three main disadvantages when com- 
pared with P95. Firstly, it has been taken in 4-CCD mode. 
The energy resolution of the SISs in 2-4 CCD modes has 
degraded by a factor of four with time with respect to 
1-CCD mode observations due to the increase of RDD^". 
RDD effects were already not negligible for 4-CCD mode 
observations at the time P93 was taken (Rasmussen et al. 
1994). Secondly, the effective exposure time is relatively 
short, around four times smaller than for P95; Thirdly, 
the cluster's central regions are mainly extended over two 
chips in each spectrometer (SOCl & S0C2 and S1C3 & 
SICO). These two chips have different overall gain off-sets, 
requiring a more elaborated correction for gain variations, 
thus increasing the uncertainties of the measured redshifts. 
The first two of the above mentioned problems cause, ef- 
fectively, larger statistical uncertainties when determining 
velocities. The last problem is more serious. If not taken 
into account, interchip gain difference can be much larger 
than the velocities one is trying to measure, thus being 
able to artificially create, erase, or amplify velocity gradi- 
ents. 

We corrected for interchip gain variations as follows: 
since the cluster's center is located very near the sepa- 
ration between chips, we selected a region with 3' radius 
centered on the X-ray peak and analyzed the spectra of 
the two parts of that region that fell on each separate 
chip; 3' at Centaurus distance corresponds to ~57 h^Q kpc. 
At that distance from the cluster center we assume that 
any velocity gradient that we measure between the two 
semi-regions is purely due to interchip gain variations. If 
there were true velocity differences detectable by the spec- 
trometers in that region (^600 km s^^) the gas would be 
gravitationally unbound. By comparing the velocity dif- 
ference between these two regions we found the correction 
factors that should be applied to each CCD so that the 
gain difference between chips did not affect the velocity 
measurements. To join data from both spectrometers (to 
improve statistics) we corrected for the global gain differ- 
ences for each pair of chips in each spectrometer. In order 
to bring all redshift values from all four chips to a unique 
solution within the central region we determined the gain 
correction factors that should be applied simultaneously to 
SOCl, S0C2, S1C3 & SICO as 1.0%, 1.2%, 0.6% & 1.0%, 
respectively. 

After applying the gain correction described above we 
selected extraction regions analogous to those used for 
P95, i.e. with 3' radius at 5' from the center with a ~ 
20° rotation to compensate for the roll angle difference 
between the two pointings. If one of these regions crossed 
the chip boundaries and therefore could be separated in 
two parts we ignored the part that had the smaller num- 
ber of counts (the improvement in photon statistics does 
not justify the uncertainties generated by the gain correc- 
tion procedures). 

The results are shown if Figure 6 and listed in Table 
5. The temperature and abundance distributions are con- 
sistent with, though less significant than, those obtained 

heasarc.gsfc.nasa.gov/docs/asca/newsletters/sis_calibration5.html 
heasarc.gsfc.nasa.gov/docs/asca/ newsletters/sis_performance2.html 

Within the extraction regions considered in the analysis of the central regions (circles of 3' radius 5' from the center) some photons fell in 
the top two chips (in detector coordinates) of each spectrometer (SOCO & S0C3 and S1C2 & SlCl). However, these regions are small and the 
gain correction procedure described above would create much higher uncertainties than what we would gain statistically by including those 
regions. Therefore, we consider photons detected only within the two chips described in the previous paragraph. 



for P95. The azimuthal temperature gradient outside the 
cooling flow region is clearly seen. The maximum temper- 
ature is reached towards the Cen 45 group (T=3.49±0.25 
keV), the temperature within the central 3' is 2.54tQ'Q4 
keV and the central abundance gradient is marginally vis- 
ible. The large errors associated with velocity measure- 
ments in P93 are due to: 1) the effective exposure being 
~ 4 times smaller than in P95; 2) the observations being 
taken in 4-CCD mode, which worsens the energy resolu- 
tion (^3.5% at 6.0 keV); 3) the correcting procedure for 
interchip gain variations involves assumptions that may 
be oversimplified (e.g. the non-existence of true gas veloc- 
ity differences within a 3' radius. Nevertheless, a velocity 
gradient is also seen roughly in the same orientation as 
that detected in P95! In particular, P3 & P4 have signif- 
icantly (> 68% confidence) lower redshifts than P6. The 
velocity difference correspondent to regions P4 & P6 is 
AV'^ (1.8 ±1.5) X lO'^fcm s~^ (1-cr errors), consistent with 
the velocity measurements from P95. 



6. DISCUSSION 

Spatially resolved spectral analysis of ASCA SIS & 1 
data carried out in this work strongly indicates the pres- 
ence of significant azimuthal temperature, metal abun- 
dance and velocity gradients in the ICM of the Centaurus 
cluster. The temperatiire and metal abundance profiles 
derived in this work indic;ate the presence of a very mild 
cooling flow and a central abundance enhancement, con- 
sistent with those found in previous works. Both temper- 
ature and abundance gradients are steeper (with different 
signs) towards the Southeast direction, where Cen 45 is 
located. The increase in temperature and mild dilution of 
metal abundances towards Cen 45 is consistent with the 
hypothesis that Cen 45 is in the process of merging with 
the main body of the Centaurus cluster (Cen 30). 

The velocity distribution shows a significant gradient 
along a direction roughly perpendicular to the direction of 
Cen 45 subgroup. The amplitude of the velocity gradi- 
ent (along the P4-P6 direction) can be as high as ~ 25 hso 
km s^^kpc^^. The velocity gradient is significantly (>90% 
confidence) detected with both SISs on-board ASCA. The 
velocity gradient is also seen in a different ASCA obser- 
vation performed 2 years earlier, although with smaller 
significance. We have shown in the previous sections that 
the velocity gradient is not caused by gain fluctuations 
across the SIS chips, thus suggesting the presence of large- 
scale bulk motions of the intracluster gas in this cluster. 

A tentative global velocity distribution for the central 
regions of the Centaurus cluster that summarizes the re- 
sults obtained in this paper is plotted in Figure 7 and listed 
in Table 6. The data points show a weighted average ve- 
locity departure from the median AV for all cases that we 
have analyzed previously, i.e.. SIS & 1 for P95, SIS & 
1 for P95 without the FeK line complex and SIS & 1 for 
P93. The errors include the fltting errors and a constant 
gain uncertainty corresponding to the I-ct intrachip gain 
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the two sub-systems (Ricker 1998), and thirdly, galaxies 
in the two sub-systems in Centaurus are still in a clearly 
segregated state. 

Even though it has been suggested that mergers erase 
cooling flows (e.g. Edge, Steward & Fabian 1992, Roct- 
tiger et al. 1993) and central abundance gradients, recent 
simulations of cooling-flow-cluster mergers show that, even 
in head-on mergers, the bow shock that appears in be- 
tween the colliding systems can protect the gas from mix- 
ing and cooling flows c;an survive mergers depending on 
the produced ram-pressure of the gas in the infalling clus- 
ter (Gomez et al. 2000, see also Fabian & Daines 1991). 
Cooling flow survival in mergers is predicted to be even 
more probable in off-center mergers (Gomez et al. 2000). 
Further, large scale merger-induced gas rotation can last 
for a few crossing times (Roettiger & Flores 2000), which 
suggests a more likely scenario to explain the gas velocity 
distribution in Centaurus, that is, the velocity gradient in 
Centaurus pre-existed the current merger and is not con- 
nected to the recent infall of Cen 45. In this scenario the 
observed velocity gradient was generated by some previ- 
ous off-center merger, that happened > 0.2-0.6 h^g^ Gyr 
(circulation time) ago and did not destroy either the cool- 
ing flow or the central abundance gradient This past 
merger event could have been a reasonably isolated event, 
in which case the temperature asymmetry in the maps of 
Churazov et al.(1999) could actually be due to the effects 
of a current head-on collision with a prc;viously rotating 
cluster. If this is the case, Centaurus may provide a crucial 
observational basis for future comparison to (currently un- 
available) numerical simulations of cluster mergers where 
the components have a previously imparted intrinsic an- 
gular momentum. 

However, the velocity structure observed could be also 
created by a past series of matter accretion events from 
fllamentary structures, in the intersection of which clus- 
ters are believed to form. The poor spatial resolution of 
the spectrometers on-board ASCA docs not allow us to 
construct a detailed velocity map that would discriminate 
between a more or less systematic ICM "rotation" from 
a more chaotic ICM "stormy weather" as predicted by 
Burns et al. 1999. Multiple off-center long exposure ve- 
locity measurements of intracluster gas with the Chandra. 
and XMM-Newton satellites will allow us to reduce the 
systematic uncertainties related to gain variations, and 
conscqiiently, to determine ICM velocities in the central 
and outer regions more precisely. This, combined with 
more speciflc cluster-cluster merger simulations will allow 
us to strongly constrain the current evolutionary stage of 
the Centaurus cluster. 

7. SUMMARY 

We have used, for the first time, spatially resolved spec- 
troscopy for determination of ICM radial- velocities in the 
Centaurus cluster. The main results of our velocity anal- 
ysis are: 

1. The gas velocity distribution at > 8' determined 

^•^ It should be noted that bulk velocities as high as these can be, in principle, measured also through high resolution differential measurements 
of the Sunyaev-Zel'dovich effect across the line of maximum velocity difference, especially in cool clusters such as Centaurus, where the thermal 
component is weaJjly dominant below 300 GHz and {^-)cmbr ~ lO"** — 10~^. 

After this paper was submitted we learned that an extensive gas temperature analysis of multiple pointings of the Centaurus cluster by 
Purusho et al. (2001) also supports this scenario. 



fluctuation described in the previous sections, i.e., 0.28%. 
The direction towards Cen 45 is shown by the vertical lines 
and the uncertainties of the velocity measurements for the 
very central regions is indicated by the horizontal lines. 
For illustration, we fit a simple solid body rotation fimc- 
tion to the data and obtain a best-fit maximum rotational 
velocity of 1.59±0.32 x 10^ km s'^ 

Similarly to the case of the Perseus cluster, the large 
scale X-ray elongation is near perpendicular to the direc- 
tion where the velocity gradient is maximum. However, in 
the case of Centaurus, the elongation towards Cen 45 can 
be directly associated with the elongation of the gravita- 
tional well along the collision axis. At smaller scales (^ 
10') there is evidence for isophotal twisting (Mohr, Fab- 
ricant & Geller 1993) and the major axis of the surface 
brightness elongation shifts ~ 50° towards the direction of 
maximum gas velocity gradient. 

One way of generating such large angular momentum in 
clusters is through cluster-cluster mergers. Several recent 
numerical simulations of off-center cluster mergers predict 
residual intracluster gas bulk rotation with velocities of a 
few thousand km s^^ (Ricker 1998; Roettiger, Stone, & 
Mushotzky 1998; Takizawa 1999, 2000; Roettiger & Flo- 
res 2000) indicating that a high fraction of the merging en- 
ergy can be transferred to gas rotation (Pearce, Thomas, & 
Couchman, 1994). In the case of Centaurus there is grow- 
ing evidence indicating that this cluster is in the process 
of merging with the galaxy group Cen 45 including: (1) 
the galaxy radial velocity bimodality detected by Lucey 
et al. (1986a,b), (2) the large scale gas radial velocity bi- 
modality shown in this work (CIS analysis); (3) the large 
scale temperature gradient detected towards NGC 4709, 
which is expected from heating due to interaction between 
the two subgroups (Churazov et al. 1999); (4) the small 
scale azimuthal temperature and abundance gradients de- 
tected towards Cen 45 (this work); (5) the morphological 
segregation both in velocity and physical space: Cen 30 is 
dominated by early type galaxies and Cen 45 by late types 
(Stein, .Icrjcn & Federspiel 1997); 

Simulations of head-on mergers show that the core of the 
primary cluster changes from being spherical before core- 
crossing to an elongated shape along the merger direction 
after core-crossing (Gomez 2000). This is not observed 
in the core of Cen 30 (the X-ray elongation in the core 
is actually near perpendicular to the collision axis). In a 
head-on collision with a smaller mass sub-unit, the core of 
the main cluster may display cooler regions towards the 
side opposite to the strongest bow shock (at the boundary 
of the two systems) shortly (0.25 Gyr) after core crossing 
(e.g. Gomez 2000). This would be consistent with the 
smooth azimuthal temperature gradient along the direc- 
tion towards Cen 45 (PI, PO, P5). However, this scenario, 
i.e., we are in a post-core crossing epoch, faces several diffi- 
culties. Firstly, this is a very short period of time in order 
to establish full rotation. Secondly, even if rotation is not 
fully established significant gas bulk motions are seen only 
after the boundary shock has encompassed the cores of 
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with the GISs shows that the region associated 
with Cen 45 has a radial velocity higher than 
that of the rest of the cluster by ~ 1.77±0.96 x 
10^ km s~^. This is consistent with the optically 
determined velocity differences between Cen 45 
and Cen 30 of ~ 1.35±0.19 x 10^ km s'^ (Stein et 
al. 1997). 

2. We found azimuthal temperature and abundance 

gradients in the cluster's central regions (< 8'). 
The azimuthal temperature variations are in 
agreement with the temperature maps determined 
previously by other authors (Churazov et al. 
1999). Both temperature and metal abundance 
variations seem to be directly related to the infall 
of the Cen 45 subgroup. 

3. The gas velocity distribution at < 8' determined 
with the SISs reveals a significant velocity 
gradient (>99.9% confidence) along a direction 
ncar-pcrpcndicular to the direction towards Ccn45. 
This detection is found to be significant in both 
SIS & 1 separately and also in another ASCA 
pointing taken two years earlier. The velocity 
gradient does not show any clear correlation with 



the temperature or metal abundance gradients. 
It is not due to spectral modeling effects or to 
intrachip gain fluctuations. The velocity gradient 
is consistent with gas bulk rotation with circular 
velocity 1.59±0.32 x 10"^ km s~^, implying that 
gas kinetic energy is comparable to the thermal 
energy in the intracluster gas of Centaurus. 

4. Comparison to cluster-merger numerical simula- 
tions, available in the literature, suggests that 
the velocity gradients that we detected in the 
central regions of Centaurus are not associated 
with the current infall of the Cen 45 subgroup. 
Therefore, we believe that it is more likely that the 
central velocity gradient is remnant of one or more 
previous off-center merger events. 

We would like to thank J. Irwin for the many helpful 
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useful suggestions. We particularly thank K. Mukai for 
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FIGURE CAPTIONS 

Fig. 1. — Distribution of the spatial regions for the pointings analyzed in this work. Regions labeled with a Pj, where i=0-8, denote the 
set of extraction regions used in the SIS &; 1 analysis of the central regions for P95 and P93. The radius of each circular region is 3' and 
they are 5' away from the cluster's center. The larger regions denoted J12, J34, J56 and J78 are combinations two-by-two of circular regions 
12' away with 5' radius used for the GIS velocity analysis of the outer regions of Centaurus. For the same analysis a circular region 11' away 
with 8' radius was also selected and is denoted by C45. The approximate position of NGC 4709 is towards the bottom left of C45. Surface 
brightness contours of Centaurus arc overlaid. 

Fig. 2. — a)Azimuthal distribution of temperature (TOP), metal abundance (MIDDLE), and redshift (BOTTOM) as a function of the 
position angle (radians) obtained from the analysis of SIS data for P95. The first data point from the left for all plots corresponds to PI, 
increasing to P8 (last). For all plots dark solid horizontal lines represent the 90% confidence limits for the central region (PO). The region 
between vertical lines indicates the approximate direction of Cen 45. Errors for all plots are 90% confidence. A velocity difference equivalent 
to 2700 km s~^ is also indicated by horizontal dashed lines. b)Same as (a) but for SIS 1 

Fig. 3. — Confidence contour plot for two interesting parameters (redshifts) for regions P4 & P6 using SIS & 1 data for P95. The two 
contours correspond to 90% and 99% confidence levels. The line of equal redshifts is also indicated. The contours are found from simultaneous 
spectral fittings of four data groups (P4 SIS O&l and P6 SIS O&l). The redshifts parameters for both instruments axe tied together for the 

spectral fittings of the same region. 

Fig. 4. — Probability of finding redshift differences equal or greater than what we observed for the real data (P4xP6) by chance as a 
function of the standard deviation of intrachip gain variations for the SISs. The cTgain for SIS is assumed to be the same as that for SIS 1 
in this plot for illustration purposes. The actual value of again used in this paper (^0.28%) is also indicated. 

Fig. 5. a)Azimuthal distribution of redshifts as a function of the position angle (radians) obtained from the spectral fittings of SIS data 
for P95 when the FeK line complex energy range was excluded. First data point from the left for all plots corresponds to PI, increasing to 
P8 (last). For all plots dark solid horizontal lines represent the 90% confidence limits for the central region (PO). The region between vertical 
lines indicate the approximate direction of Cen 45. Errors for all plots are 90% confidence. A velocity difference equivalent to 6300 km is 
also indicated by horizontal dashed lines. b)Same as (a) but for SIS 1 and the horizontal dashed lines indicated a velocity difference equivalent 
to 4300 km s-i. 

Fig. 6. — Azimuthal distribution of Temperature (TOP), Metal Abundance (MIDDLE), and Redshift (BOTTOM) as a function of the 
position angle (radians) obtained from the simultaneous spectral fittings of SIS O&l data for P93. The first data point from the left for all 
plots corresponds to PI, increasing to P8 (last). The errorbars indicated show 90% confidence errors for temperature and abundances and 
68%(1(t) for redshifts. The dark solid horizontal lines represent the confidence limits for the central region (PO). The region between vertical 
lines indicate the approximate direction of Cen 45. A velocity difference equivalent to 1800 km is also indicated by horizontal dashed 
lines. 

Fig. 7. — Average aaimuthal gas radial velocity distribution for the central regions of the Centaurus cluster. Data points indicate the 
weighted average of the departure from the median. The errors also include a la gain variation of 0.28% (typical for the SIS chips at 1995) 
for all cases (all instruments, all pointings) analyzed in this work. The first data point from the left for all plots corresponds to PI, increasing 
to P8 (last). The dark solid horizontal lines represent the 68% confidence limits for the central region (PO). The region between the two 
vertical lines indicates the approximate direction of Cen 45. The solid line indicate the bestfit for a simple solid-body cosine function with a 
corresponding maximum rotational velocity of ~ 1600±300km s~^. 



Table 1 
X-RAY Observations 



Pointing 


Sequence 


Date 


RA DEC 


CIS EXP*^ 


SIS EXPb 




Number 


Observed 


(2000) (2000) 


(ksec) 


(ksec) 


P95 


83026000 


1995-07-19 


12h48m21.60s -41°19'44.8" 


56.4 


67.9 


P93 


80032000 


1993-06-30 


12h48m24.00s -41°17'16.4" 


14.7 


16.6 



"^Effective Exposure (Average for GIS 2 & 3) 
''Effective Exposure (Average for SIS & 1) 



Table 2 

GIS Spectral Fittings for P95''' 



Region 


Temperature'^ 
(keV) 


Abund 
(Solar) 


Redshift 
(10-2) 




C45 

A11-C45 

JO 

J12 

J34 

J56 

J78 


r 47+0.28 
^•+' -0.27 
o 44+0.03 
"^•^+-0.03 
9 70+0.02 
^- '"^-0.02 
(^7+0.09 
•J-" ' -0.08 

4 16+°-i3 

^■J^"-0.13 

^•0' -0.17 
oq+0.10 


fin+" "^ 

'J-OU-o.OO 
n 77+0.03 
^- ' ' -0.03 

1 09+005 



L'-o'i-o.oe 

-0.08 

n 44+0-09 

'-'•^^-0.08 

49+0.06 


q 97+1.58 
•^•^ ' -1.13 

1 n-,+0.05 

i-oi_o.oi 

1 01 +0.05 

i-oi_o.oi 

9 qr+0.37 

^••J^-0.52 

1 qq+1-07 

^•^'-'-0.08 
^•^^-0.50 


1.062 
1.206 
1.102 
0.983 
0.969 
0.969 
0.970 



'^Errors are 90% confidence 
^0.8-10.0 keV 

'^Nh fixed at the Galactic value 



Table 3 
Spectral Fittings for P95''''''° 



Region 


T(SIS0)'i'^ 


T(SISl)'' 


Abund (SISO) 


Abund (SISl) 


z (SISO) 
(10-2) 


z (SISl) 

(10-2) 




(keV) 


(keV) 


(Solar) 


(Solar) 


PO 


9 P7+0.03 

^■"'-0.03 


9 pq+o.oe 

^•"^-0.05 


-r iq+0-06 
^^"^-0.03 


-1 9q+0.07 

■■-•^•J-o.oe 


1 q9+0.02 

-■-••J^-o.oi 




PI 


2 04+0.06 

^•^^-0.06 


9 7/1+0.08 
^- '*-0.06 


1 11+0-12 
J^'-'^-'^-O.lO 


1 10+0-15 

J^--l"-0.12 


1 00+0. 15 
^•'J'-'-O.lO 


1 ^8+0-38 
i.OO_Q 23 


P2 


q -r 1 +0.07 
"^■^^-0.06 


02+0 07 
^■^^-0.07 


96+0 09 

^■^"-0.09 


1 11+0-12 
■'■•-'-^-0.11 


-1 qq + 0.16 

J^-'J'-'-O.ll 




P3 


q /,q + 0.09 

'-'•^'J-0.08 


3.54^ro^ 


n nq+0.08 

^•550-0.08 




o.88t°;i;^ 


o.88i°;J^ 


P4 


rq+0.08 




0.89to-09 


0.84tO;iJ 


0.88tS3;S« 


n sr. +0.34 


P5 


c«+0.10 
^•^"-0.09 




0.83±o-09 


0.87t0;i3 


-I qq+0.12 

^•'-''-'-0.09 


1 00+0. 16 
^^•'-''-'-0.51 


P6 


q qK+0.09 
"^■•JO-O.OO 


q q7 + 0.11 

"-"••J' -0.11 




qq+O-14 

"•^■^-0.13 


1 70+0-10 
-■-•'■^-0.04 




P7 


9 qR+0.08 
^■^"-0.08 


9 qfj+o.io 

^■^"-0.09 


Q4+0-05 
^■^^-0.12 


Q'-i+o i'' 
'-'•^•J-0.12 


1 70+0.10 
J^-'^-0.06 


0.87tro§ 


PS 


9 qq+0.07 


9 88+0.08 
^•»»_0.07 


1 02+0- 

J^-^^-0.06 


1 m+012 


1 70+0.07 
■•■•'0-0.04 





'^Errors arc 90% confidence level 
^0. 8-10.0 keV 
^X^ ~ 1-1 

'^'Nh fixed at the Galactic value 
Global redshifts boosted up by 0.9% 
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Table 4 
Spectral Fittings for 



Region T(SISO)'i'" T(SIS1)<^ Abund (SISO) Abund (SISl) z (SISO) z (SISl) 
(keV) (keV) (Solar) (Solar) (10"^) (10"^) 





2. 




2, 


77+0.03 
'-0.07 


1, 


46+0.06 

•+"-0.06 


1, 


/iq+0.07 
■^'-'-0.08 


2, 


71 +0.02 

• ' -"--0.00 


2, 


07+0.04 

■"5' -0.07 




3, 




2, 


OR+0.05 
■^0-0.11 


1, 


1Q+0.15 
•^"-0.16 


1, 


01 +0.25 
•"^^-0.20 


2 


90+0.47 
■^"^-o.ig 


2, 


qc+0.07 
■"-"^-0.09 


P2FeL 


3, 


01 +0.05 
■^^-0.10 


3, 


1 o+O.lO 
■^°-0.09 


1, 


iq+o.i7 

■-■-^-0.14 


1, 


00+0.20 
•^"-0.18 


1, 


cq+0.69 
■">5_o.08 


2, 


qc+0.06 
■"5"-0.06 


P^FeL 


3, 


rco+0.11 
■oO_Q 12 


3, 


09+0. 07 
■°^-0.13 


1, 


■^^-0.13 


1, 


t:;t:;+0.07 
■■30-0.23 


1, 


iq+0.67 
■-^^-0.07 


2, 


09+0.09 
■"5^-0.55 


P^FeL 


3, 


70+0.13 
■'^-0.11 


4, 


^^+0.17 
■^^-0.16 


1, 


1 0+0. 16 
■^"^-0.15 


1, 


90+0.28 
■^"5-0.18 


1, 


71+0.14 
■ ' -^-0.19 


1, 


qq+0.34 
■"5"5-0.50 


P^FeL 


3, 


71+0.14 
■ ' .'--0.14 


3, 


qq+0.19 

■^^-0.18 


1, 


r,i+0.20 

•^+-0.16 


1, 


Qq+0.14 

•"■J-0.21 


2, 


r7+0.82 
•0 ' -0.40 


1, 


q7 + 0.91 

■"5' -0.14 


P^FeL 


3, 


59+0.10 

■^^-0.10 


3, 


cc+0.12 
■^^-0.19 


1, 


20+0.23 
•^■^-0.23 


1, 


9 c +0.25 
■^°-0.26 


3 


fjn+0.47 
■""-0.17 


2, 


79+0.18 
■'^-0.98 


P^FeL 


3, 


1 1+0.14 
■-^^-0.12 


3, 


r,q + 0.15 

■^•5-0.15 


1, 


1 7+0.24 
' -0.17 


1, 


04+O.2S 
■"5^-0.24 


3 


07+0.54 
■"' -0.51 


2, 


54+0.63 
■■J+-0.44 


PSFeL 


3, 


10+0.07 
•J- "5-0.09 


3, 


07+0.13 
•U' -0.11 


1. 


45+0.19 


1. 


qq+0.27 
■ "5 "5-0. 18 


2. 


-?q + 0.22 

■""5-0.08 


2. 


90+0.09 
■"'°-0.48 



'^Errors arc 90% confidence level 
'^FeK complex excluded (0.8-6 keV) 
^X^ ~ 1-1 

'^'Nh fixed at the corresponding Galactic value 
''Global redshifts boosted up by 1.2% 



Table 5 
Spectral Fittings for 



Region 


Temp^ 


Abund " 








(keV) 


(Solar) 


(10-2) 




PO 


9 r4+0.07 

^•3^-0.04 




1.02tO;°^ 


1.19 


PI 


"5-""-0.14 




^•"^-0.15 


1.125 


P2 


on+0.14 
•5-'J"-0.15 


o.9it°;}? 


1 1 c+0.44 
l^-J^O_o.06 


1.132 


P3 


q qt;+0.16 
"5-"5O-0.15 


O.Slt";- 


^•""-0.07 


0.994 


P4 


46+0.20 


0.87tS3:}^ 


75+0-27 


1.095 


P5 


4q+o.26 

'5-^-J-0.24 


75+0.16 

'J- ' •J-0.17 


Q 77+0.79 
U. ( ( _Q 21 


1.013 


P6 


q 97+0.17 
"5-^ ' -0.16 


"■ ' ^-0.09 


1 q9+0.27 
-■-•"5^-0.22 


0.995 


P7 


q iq+0.13 
"5-J-"5-0.13 




1 11+0.31 
■'-•-'■-'--0.15 


1.055 


P8 


2.80l°;fe 


0.851H! 


^^•'^"-0.06 


1.399 



'^Nh fixed at the Galactic value 

^Errors are 68% confidence level 

'^Errors arc 90% confidence level 

Slope boosting correction for interchip gain differ- 
/ 50C1 50C2 \ _ / 1.01 1.012 \ 
'^^^ \ S1C3 SICO j ~ \ 1.006 1.01 j 



Table 6 

Weighted Average Velocity Differences for All Pointings 



Region 


AV 




(103 km s-i) 


PO 


0.07±0.45 


PI 


-0.06±0.52 


P2 


-0.63±0.47 


P3 


-1.84±0.45 


P4 


-2.04±0.54 


P5 


-0.36±0.74 


P6 


1.38±0.53 


P7 


0.71±0.61 


P8 


0.46±0.49 



'^la errors include spectral 
fitting as well as intrachip gain 
variations. 



LU 



4.0 



3.5 . 



<> 

lU 2.5 . 



LU 
O 



1.3 
1.2 
.1.1 



||1.0 
5 W 0.9 



CQ 
< 



0.8 
0.7 



0.022 

0.020 

l_ 0.018 
U. 

^ 0.016 
CO 

Q 0.014 ^ 
LU 

^ 0.012 
0.010 ■ 
0.008 • 



SISO 
P95 



I 



P1 P2 



I I 



P6 

-I- 



P7 P8 



P5 



P3 P4 



I 



-1 1 2 3 4 5 

PA (rad) 



4.0 



LU 

OC 3.5 



<> 

OC 0)3.0 

LU ^ 



LU 



LU 
O 
Z 
< 
Q 



CQ 
< 



CO 
Q 
LU 
OC 



2.5 

2.0 
1.3 

1.2 

.1.1 

iSl.O 
O 

^0.9 
0.8 
0.7 

0.022 
0.020 
0.018 
0.016 
0.014 
0.012 
0.010 
0.008 



SIS1 
P95 



P1 



P2 



P6 



P7 



P8 



P5 



PA 



P3 



-1 1 2 3 4 5 

PA (rad) 




(9d} UiMBpod 



0.010 
0.009 
0.008 
0.007 
— 0.006 
■§0.005 

n 

00.004 
°" 0.003 
0.002 
0.001 
0.000 




used 

gam 



gain 



P95 — FeL 





0.040 




0.035 


H 




u- 0.030 


X 




DS 


0.025 


LU 




OC 


0.020 




0.015 




0.010 



(a) SIS 



P1 



P2 



P3 p4 





0.035 




0.030 






LL 


0.025 


X 




CO 




Q 0.020 


LU 




CC 






0.015 




0.010 




0.005 



(b) SIS 1 



P1 



P2 



P4 



P3 



P6 



P7 



P5 



P7 



P6 



P5 



1 2 
PA (rad) 



P8 



P8 . 



LU 



4.0 



3.5 . 



<> 

S 3.0 ■ 



LU 



LU 
O 



2.5 

1.1 
1.0 



< feO.9 
Z §0.8 



CQ 
< 



0.7 
0.6 



0.020 

0.018 ^ 

l_ 0.016 
LJ. 

^ 0.014 



0.006 ■ 



SIS &1 
P93 



I 




I 



1 2 
PA (rad) 



I 



^ 1000 

en 



SIS 0&1 
P95 & P93 




12 3 

PA (rad) 



